Longitudinal change in mitochondrial heteroplasmy exhibits positive selection for deleterious variants.

A common feature of human aging is the acquisition of somatic mutations, and mitochondria are particularly prone to mutation due to their inefficient DNA repair and close proximity to reactive oxygen species, leading to a state of mitochondrial DNA heteroplasmy1,2. Cross-sectional studies have demonstrated that detection of heteroplasmy increases with participant age3, a phenomenon that has been attributed to genetic drift4-7. In this first large-scale longitudinal study, we measured heteroplasmy in two prospective cohorts (combined n=1405) at two timepoints (mean time between visits, 8.6 years), demonstrating that deleterious heteroplasmies were more likely to increase in variant allele fraction (VAF). We further demonstrated that increase in VAF was associated with increased risk of overall mortality. These results challenge the claim that somatic mtDNA mutations arise mainly due to genetic drift, instead demonstrating positive selection for predicted deleterious mutations at the cellular level, despite an negative impact on overall mortality.

A new variant in the ZCCHC8 gene: diverse clinical phenotypes and expression in the lung.

Introduction: Pulmonary fibrosis is a severe disease which can be familial. A genetic cause can only be found in ∼40% of families. Searching for shared novel genetic variants may aid the discovery of new genetic causes of disease. Methods: Whole-exome sequencing was performed in 152 unrelated patients with a suspected genetic cause of pulmonary fibrosis from the St Antonius interstitial lung disease biobank. Variants of interest were selected by filtering for novel, potentially deleterious variants that were present in at least three unrelated pulmonary fibrosis patients. Results: The novel c.586G>A p.(E196K) variant in the ZCCHC8 gene was observed in three unrelated patients: two familial patients and one sporadic patient, who was later genealogically linked to one of the families. The variant was identified in nine additional relatives with pulmonary fibrosis and other telomere-related phenotypes, such as pulmonary arterial venous malformations, emphysema, myelodysplastic syndrome, acute myeloid leukaemia and dyskeratosis congenita. One family showed incomplete segregation, with absence of the variant in one pulmonary fibrosis patient who carried a PARN variant. The majority of ZCCHC8 variant carriers showed short telomeres in blood. ZCCHC8 protein was located in different lung cell types, including alveolar type 2 (AT2) pneumocytes, the culprit cells in pulmonary fibrosis. AT2 cells showed telomere shortening and increased DNA damage, which was comparable to patients with sporadic pulmonary fibrosis and those with pulmonary fibrosis carrying a telomere-related gene variant, respectively. Discussion: The ZCCHC8 c.586G>A variant confirms the involvement of ZCCHC8 in pulmonary fibrosis and short-telomere syndromes and underlines the importance of including the ZCCHC8 gene in diagnostic gene panels for these diseases.

Disruption of TUFT1, a Desmosome-Associated Protein, Causes Skin Fragility, Woolly Hair, and Palmoplantar Keratoderma.

Desmosomes are dynamic complex protein structures involved in cellular adhesion. Disruption of these structures by loss-of-function variants in desmosomal genes leads to a variety of skin- and heart-related phenotypes. In this study, we report TUFT1 as a desmosome-associated protein, implicated in epidermal integrity. In two siblings with mild skin fragility, woolly hair, and mild palmoplantar keratoderma but without a cardiac phenotype, we identified a homozygous splice-site variant in the TUFT1 gene, leading to aberrant mRNA splicing and loss of TUFT1 protein. Patients' skin and keratinocytes showed acantholysis, perinuclear retraction of intermediate filaments, and reduced mechanical stress resistance. Immunolabeling and transfection studies showed that TUFT1 is positioned within the desmosome and that its location is dependent on the presence of the desmoplakin carboxy-terminal tail. A Tuft1-knockout mouse model mimicked the patients' phenotypes. Altogether, this study reveals TUFT1 as a desmosome-associated protein, whose absence causes skin fragility, woolly hair, and palmoplantar keratoderma.

Fat metabolism is associated with telomere length in six population-based studies.

Telomeres are repetitive DNA sequences located at the end of chromosomes, which are associated to biological aging, cardiovascular disease, cancer and mortality. Lipid and fatty acid metabolism have been associated with telomere shortening. We have conducted an in-depth study investigating the association of metabolic biomarkers with telomere length (LTL). We performed an association analysis of 226 metabolic biomarkers with LTL using data from 11 775 individuals from six independent population-based cohorts (BBMRI-NL consortium). Metabolic biomarkers include lipoprotein lipids and subclasses, fatty acids, amino acids, glycolysis measures and ketone bodies. LTL was measured by quantitative polymerase chain reaction or FlowFISH. Linear regression analysis was performed adjusting for age, sex, lipid-lowering medication and cohort-specific covariates (model 1) and additionally for body mass index (BMI) and smoking (model 2), followed by inverse variance-weighted meta-analyses (significance threshold Pmeta = 6.5 × 10-4). We identified four metabolic biomarkers positively associated with LTL, including two cholesterol to lipid ratios in small VLDL (S-VLDL-C % and S-VLDL-CE %) and two omega-6 fatty acid ratios (FAw6/FA and LA/FA). After additionally adjusting for BMI and smoking, these metabolic biomarkers remained associated with LTL with similar effect estimates. In addition, cholesterol esters in very small VLDL (XS-VLDL-CE) became significantly associated with LTL (P = 3.6 × 10-4). We replicated the association of FAw6/FA with LTL in an independent dataset of 7845 individuals (P = 1.9 × 10-4). To conclude, we identified multiple metabolic biomarkers involved in lipid and fatty acid metabolism that may be involved in LTL biology. Longitudinal studies are needed to exclude reversed causation.

Impact of SNP microarray analysis of compromised DNA on kinship classification success in the context of investigative genetic genealogy.

Single nucleotide polymorphism (SNP) data generated with microarray technologies have been used to solve murder cases via investigative leads obtained from identifying relatives of the unknown perpetrator included in accessible genomic databases, an approach referred to as investigative genetic genealogy (IGG). However, SNP microarrays were developed for relatively high input DNA quantity and quality, while DNA typically obtainable from crime scene stains is of low DNA quantity and quality, and SNP microarray data obtained from compromised DNA are largely missing. By applying the Illumina Global Screening Array (GSA) to 264 DNA samples with systematically altered quantity and quality, we empirically tested the impact of SNP microarray analysis of compromised DNA on kinship classification success, as relevant in IGG. Reference data from manufacturer-recommended input DNA quality and quantity were used to estimate genotype accuracy in the compromised DNA samples and for simulating data of different degree relatives. Although stepwise decrease of input DNA amount from 200 ng to 6.25 pg led to decreased SNP call rates and increased genotyping errors, kinship classification success did not decrease down to 250 pg for siblings and 1st cousins, 1 ng for 2nd cousins, while at 25 pg and below kinship classification success was zero. Stepwise decrease of input DNA quality via increased DNA fragmentation resulted in the decrease of genotyping accuracy as well as kinship classification success, which went down to zero at the average DNA fragment size of 150 base pairs. Combining decreased DNA quantity and quality in mock casework and skeletal samples further highlighted possibilities and limitations. Overall, GSA analysis achieved maximal kinship classification success from 800 to 200 times lower input DNA quantities than manufacturer-recommended, although DNA quality plays a key role too, while compromised DNA produced false negative kinship classifications rather than false positive ones.

Somatic TARDBP variants as a cause of semantic dementia.

The aetiology of late-onset neurodegenerative diseases is largely unknown. Here we investigated whether de novo somatic variants for semantic dementia can be detected, thereby arguing for a more general role of somatic variants in neurodegenerative disease. Semantic dementia is characterized by a non-familial occurrence, early onset (<65 years), focal temporal atrophy and TDP-43 pathology. To test whether somatic variants in neural progenitor cells during brain development might lead to semantic dementia, we compared deep exome sequencing data of DNA derived from brain and blood of 16 semantic dementia cases. Somatic variants observed in brain tissue and absent in blood were validated using amplicon sequencing and digital PCR. We identified two variants in exon one of the TARDBP gene (L41F and R42H) at low level (1-3%) in cortical regions and in dentate gyrus in two semantic dementia brains, respectively. The pathogenicity of both variants is supported by demonstrating impaired splicing regulation of TDP-43 and by altered subcellular localization of the mutant TDP-43 protein. These findings indicate that somatic variants may cause semantic dementia as a non-hereditary neurodegenerative disease, which might be exemplary for other late-onset neurodegenerative disorders.

Short-Term, Combined Fasting and Exercise Improves Body Composition in Healthy Males.

Fasting enhances the beneficial metabolic outcomes of exercise; however, it is unknown whether body composition is favorably modified on the short term. A baseline-follow-up study was carried out to assess the effect of an established protocol involving short-term combined exercise with fasting on body composition. One hundred seven recreationally exercising males underwent a 10-day intervention across 15 fitness centers in the Netherlands involving a 3-day gradual decrease of food intake, a 3-day period with extremely low caloric intake, and a gradual 4-day increase to initial caloric intake, with daily 30-min submaximal cycling. Using dual-energy X-ray absorptiometry analysis, all subjects substantially lost total body mass (-3.9 ± 1.9 kg; p < .001) and fat mass (-3.3 ± 1.3 kg; p < .001). Average lean mass was lost (-0.6 ± 1.5 kg; p < .001), but lean mass as a percentage of total body mass was not reduced. The authors observed a loss of -3.9 ± 1.9% android fat over total fat mass (p < .001), a loss of -2.2 ± 1.9% gynoid over total fat mass (p < .001), and reduced android/gynoid ratios (-0.05 ± 0.1; p < .001). Analyzing 15 preselected single-nucleotide polymorphisms in 13 metabolism-related genes revealed trending associations for thyroid state-related single-nucleotide polymorphisms rs225014 (deiodinase 2) and rs35767 (insulin-like growth factor1), and rs1053049 (PPARD). In conclusion, a short period of combined fasting and exercise leads to a substantial loss of body and fat mass without a loss of lean mass as a percentage of total mass.

Reduced penetrance of pathogenic ACMG variants in a deeply phenotyped cohort study and evaluation of ClinVar classification over time.

PURPOSE: We studied the penetrance of pathogenically classified variants in an elderly Dutch population from the Rotterdam Study, for which deep phenotyping is available. We screened the 59 actionable genes for which reporting of known pathogenic variants was recommended by the American College of Medical Genetics and Genomics (ACMG), and demonstrate that determining what constitutes a known pathogenic variant can be quite challenging. METHODS: We defined "known pathogenic" as classified pathogenic by both ClinVar and the Human Gene Mutation Database (HGMD). In 2628 individuals, we performed exome sequencing and identified known pathogenic variants. We investigated the clinical records of carriers and evaluated clinical events during 25 years of follow-up for evidence of variant pathogenicity. RESULTS: Of 3815 variants detected in the 59 ACMG genes, 17 variants were considered known pathogenic. For 14/17 variants the ClinVar classification had changed over time. Of 24 confirmed carriers of these variants, we observed at least one clinical event possibly caused by the variant in only three participants (13%). CONCLUSION: We show that the definition of "known pathogenic" is often unclear and should be approached carefully. Additionally variants marked as known pathogenic do not always have clinical impact on their carriers. Definition and classification of true (individual) expected pathogenic impact should be defined carefully.

Exome Sequencing Analysis Identifies Rare Variants in ATM and RPL8 That Are Associated With Shorter Telomere Length.

Telomeres are important for maintaining genomic stability. Telomere length has been associated with aging, disease, and mortality and is highly heritable (∼82%). In this study, we aimed to identify rare genetic variants associated with telomere length using whole-exome sequence data. We studied 1,303 participants of the Erasmus Rucphen Family (ERF) study, 1,259 of the Rotterdam Study (RS), and 674 of the British Heart Foundation Family Heart Study (BHF-FHS). We conducted two analyses, first we analyzed the family-based ERF study and used the RS and BHF-FHS for replication. Second, we combined the summary data of the three studies in a meta-analysis. Telomere length was measured by quantitative polymerase chain reaction in blood. We identified nine rare variants significantly associated with telomere length (p-value < 1.42 × 10-7, minor allele frequency of 0.2-0.5%) in the ERF study. Eight of these variants (in C11orf65, ACAT1, NPAT, ATM, KDELC2, and EXPH5) were located on chromosome 11q22.3 that contains ATM, a gene involved in telomere maintenance. Although we were unable to replicate the variants in the RS and BHF-FHS (p-value ≥ 0.21), segregation analysis showed that all variants segregate with shorter telomere length in a family. In the meta-analysis of all studies, a nominally significant association with LTL was observed with a rare variant in RPL8 (p-value = 1.48 × 10-6), which has previously been associated with age. Additionally, a novel rare variant in the known RTEL1 locus showed suggestive evidence for association (p-value = 1.18 × 10-4) with LTL. To conclude, we identified novel rare variants associated with telomere length. Larger samples size are needed to confirm these findings and to identify additional variants.

Genome-wide Association Analysis in Humans Links Nucleotide Metabolism to Leukocyte Telomere Length.

Leukocyte telomere length (LTL) is a heritable biomarker of genomic aging. In this study, we perform a genome-wide meta-analysis of LTL by pooling densely genotyped and imputed association results across large-scale European-descent studies including up to 78,592 individuals. We identify 49 genomic regions at a false dicovery rate (FDR) < 0.05 threshold and prioritize genes at 31, with five highlighting nucleotide metabolism as an important regulator of LTL. We report six genome-wide significant loci in or near SENP7, MOB1B, CARMIL1, PRRC2A, TERF2, and RFWD3, and our results support recently identified PARP1, POT1, ATM, and MPHOSPH6 loci. Phenome-wide analyses in >350,000 UK Biobank participants suggest that genetically shorter telomere length increases the risk of hypothyroidism and decreases the risk of thyroid cancer, lymphoma, and a range of proliferative conditions. Our results replicate previously reported associations with increased risk of coronary artery disease and lower risk for multiple cancer types. Our findings substantially expand current knowledge on genes that regulate LTL and their impact on human health and disease.

Telomere Length and the Risk of Alzheimer's Disease: The Rotterdam Study.

There is a wide interest in biomarkers that capture the burden of detrimental factors as these accumulate with the passage of time, i.e., increasing age. Telomere length has received considerable attention as such a marker, because it is easily quantified and it may aid in disentangling the etiology of dementia or serve as predictive marker. We determined the association of telomere length with risk of Alzheimer's disease and all-cause dementia in a population-based setting. Within the Rotterdam Study, we performed quantitative PCR to measure mean leukocyte telomere length in blood. We determined the association of telomere length with risk of Alzheimer's disease until 2016, using Cox regression models. Of 1,961 participants (mean age 71.4±9.3 years, 57.1% women) with a median follow-up of 8.3 years, 237 individuals were diagnosed with Alzheimer's disease. We found a U-shaped association between telomere length and risk of Alzheimer's disease: compared to the middle tertile the adjusted hazard ratio was 1.59 (95% confidence interval (CI), 1.13-2.23) for the lowest tertile and 1.47 (1.03-2.10) for the highest tertile. Results were similarly U-shaped but slightly attenuated for all-cause dementia. In conclusion, shorter and longer telomere length are both associated with an increased risk of Alzheimer's disease in the general population.

Intestinal microbiome composition and its relation to joint pain and inflammation.

Macrophage-mediated inflammation is thought to have a causal role in osteoarthritis-related pain and severity, and has been suggested to be triggered by endotoxins produced by the gastrointestinal microbiome. Here we investigate the relationship between joint pain and the gastrointestinal microbiome composition, and osteoarthritis-related knee pain in the Rotterdam Study; a large population based cohort study. We show that abundance of Streptococcus species is associated with increased knee pain, which we validate by absolute quantification of Streptococcus species. In addition, we replicate these results in 867 Caucasian adults of the Lifelines-DEEP study. Finally we show evidence that this association is driven by local inflammation in the knee joint. Our results indicate the microbiome is a possible therapeutic target for osteoarthritis-related knee pain.

Population-specific genetic variation in large sequencing data sets: why more data is still better.

We have generated a next-generation whole-exome sequencing data set of 2628 participants of the population-based Rotterdam Study cohort, comprising 669 737 single-nucleotide variants and 24 019 short insertions and deletions. Because of broad and deep longitudinal phenotyping of the Rotterdam Study, this data set permits extensive interpretation of genetic variants on a range of clinically relevant outcomes, and is accessible as a control data set. We show that next-generation sequencing data sets yield a large degree of population-specific variants, which are not captured by other available large sequencing efforts, being ExAC, ESP, 1000G, UK10K, GoNL and DECODE.

The transcriptional landscape of age in human peripheral blood.

Disease incidences increase with age, but the molecular characteristics of ageing that lead to increased disease susceptibility remain inadequately understood. Here we perform a whole-blood gene expression meta-analysis in 14,983 individuals of European ancestry (including replication) and identify 1,497 genes that are differentially expressed with chronological age. The age-associated genes do not harbor more age-associated CpG-methylation sites than other genes, but are instead enriched for the presence of potentially functional CpG-methylation sites in enhancer and insulator regions that associate with both chronological age and gene expression levels. We further used the gene expression profiles to calculate the 'transcriptomic age' of an individual, and show that differences between transcriptomic age and chronological age are associated with biological features linked to ageing, such as blood pressure, cholesterol levels, fasting glucose, and body mass index. The transcriptomic prediction model adds biological relevance and complements existing epigenetic prediction models, and can be used by others to calculate transcriptomic age in external cohorts.

Vitamin D Receptor Polymorphisms Are Associated with Reduced Esophageal Vitamin D Receptor Expression and Reduced Esophageal Adenocarcinoma Risk.

Epidemiological studies indicate that vitamin D exerts a protective effect on the development of various solid cancers. However, concerns have been raised regarding the potential deleterious role of high vitamin D levels in the development of esophageal adenocarcinoma (EAC). This study investigated genetic variation in the vitamin D receptor (VDR) in relation to its expression and risk of Barrett esophagus (BE) and EAC. VDR gene regulation was investigated by immunohistochemistry, reverse transcriptase-polymerase chain reaction (RT-PCR) and gel shift assays. Fifteen haplotype tagging single-nucleotide polymorphisms (SNPs) of the VDR gene were analyzed in 858 patients with reflux esophagitis (RE), BE or EAC and 202 healthy controls. VDR mRNA expression was higher in BE compared with squamous epithelium. VDR protein was located in the nucleus in BE. An rs1989969T/rs2238135G haplotype was identified in the 5' regulatory region of the VDR gene. It was associated with an approximately two-fold reduced risk of RE, BE and EAC. Analysis of a replication cohort was done for BE that confirmed this. The rs1989969T allele causes a GATA-1 transcription factor binding site to appear. The signaling of GATA-1, which is regarded as a negative transcriptional regulator, could explain the findings for rs1989969. The rs2238135G allele was associated with a significantly reduced VDR expression in BE; for the rs1989969T allele, a trend in reduced VDR expression was observed. We identified a VDR haplotype associated with reduced esophageal VDR expression and a reduced incidence of RE, BE and EAC. This VDR haplotype could be useful in identifying individuals who benefit most from vitamin D chemoprevention.

Characteristics of de novo structural changes in the human genome.

Small insertions and deletions (indels) and large structural variations (SVs) are major contributors to human genetic diversity and disease. However, mutation rates and characteristics of de novo indels and SVs in the general population have remained largely unexplored. We report 332 validated de novo structural changes identified in whole genomes of 250 families, including complex indels, retrotransposon insertions, and interchromosomal events. These data indicate a mutation rate of 2.94 indels (1-20 bp) and 0.16 SVs (>20 bp) per generation. De novo structural changes affect on average 4.1 kbp of genomic sequence and 29 coding bases per generation, which is 91 and 52 times more nucleotides than de novo substitutions, respectively. This contrasts with the equal genomic footprint of inherited SVs and substitutions. An excess of structural changes originated on paternal haplotypes. Additionally, we observed a nonuniform distribution of de novo SVs across offspring. These results reveal the importance of different mutational mechanisms to changes in human genome structure across generations.

A genome-wide copy number association study of osteoporotic fractures points to the 6p25.1 locus.

BACKGROUND: Osteoporosis is a systemic skeletal disease characterised by reduced bone mineral density and increased susceptibility to fracture; these traits are highly heritable. Both common and rare copy number variants (CNVs) potentially affect the function of genes and may influence disease risk. AIM: To identify CNVs associated with osteoporotic bone fracture risk. METHOD: We performed a genome-wide CNV association study in 5178 individuals from a prospective cohort in the Netherlands, including 809 osteoporotic fracture cases, and performed in silico lookups and de novo genotyping to replicate in several independent studies. RESULTS: A rare (population prevalence 0.14%, 95% CI 0.03% to 0.24%) 210 kb deletion located on chromosome 6p25.1 was associated with the risk of fracture (OR 32.58, 95% CI 3.95 to 1488.89; p = 8.69 × 10(-5)). We performed an in silico meta-analysis in four studies with CNV microarray data and the association with fracture risk was replicated (OR 3.11, 95% CI 1.01 to 8.22; p = 0.02). The prevalence of this deletion showed geographic diversity, being absent in additional samples from Australia, Canada, Poland, Iceland, Denmark, and Sweden, but present in the Netherlands (0.34%), Spain (0.33%), USA (0.23%), England (0.15%), Scotland (0.10%), and Ireland (0.06%), with insufficient evidence for association with fracture risk. CONCLUSIONS: These results suggest that deletions in the 6p25.1 locus may predispose to higher risk of fracture in a subset of populations of European origin; larger and geographically restricted studies will be needed to confirm this regional association. This is a first step towards the evaluation of the role of rare CNVs in osteoporosis.

Recognition of scared faces and the serotonin transporter gene in young children: the Generation R Study.

BACKGROUND: Previous research highlights the significance of a functional polymorphism located in the promoter region (5-HTTLPR) of the serotonin transporter gene in emotional behaviour. This study examined the effect of the 5-HTTLPR polymorphism on emotion processing in a large number of healthy preschoolers. METHODS: The 5-HTTLPR genotype was classified in 605 children as homozygous for the short allele (SS), homozygous for the long allele (LL), or heterozygous (LS). Emotion-processing was assessed using age-appropriate computer tasks where children matched happy, sad, angry, and fearful facial expressions preceded by a shape-matching task to assess basic matching ability. RESULTS: We found that young children could differentiate between emotion categories (F = 12.1, p < .001). The effect of 5-HTTLPR genotype depended on the emotion category presented (F = 2.3, p = .031). This effect was explained by the finding that SS children were less accurate at recognising fearful faces than LL or LS children (F = 5.3, p = .005). We did not find any significant differences as a result of 5-HTTLPR genotype for happy, sad or angry expressions (p > .05). CONCLUSIONS: Results indicate that 5-HTTLPR allele status selectively impacts the processing of fearful but not other facial expressions. This pattern is already apparent in very young typically developing children. Results may signal an early vulnerability for affective problems before disorders emerge.

Serotonin transporter polymorphism moderates effects of prenatal maternal anxiety on infant negative emotionality.

BACKGROUND: Consistent with the fetal programming hypothesis, effects of maternal prenatal anxiety have been found to predict various measures of infant temperament in the early postnatal period. In recent years, a polymorphism in the serotonin transporter gene (5-HTTLPR) emerged as a moderator of diverse environmental influences on different outcomes, with individuals carrying the short allele being generally more vulnerable to adversity. METHODS: We tested whether the association between self-reported maternal anxiety at 20 weeks gestation (Brief Symptom Inventory) and mother-rated infant negative emotionality at 6 months after birth (Infant Behavior Questionnaire-Revised) would be moderated by the 5-HTTLPR in a large Dutch cohort sample (n = 1513). We hypothesized that infants carrying the 5-HTTLPR short allele would be more susceptible and therefore more affected by both low and high prenatal maternal anxiety vis-à-vis negative emotionality than other genotypes. RESULTS: Findings of a significant gene × environment interaction (B = .65, p = .01) were supportive of a vulnerability model, with infants carrying the short allele being more negatively emotional when mothers reported anxiety during pregnancy, whereas there was no difference between genotypes on negative emotionality when maternal anxiety was low. CONCLUSIONS: The association between maternal anxiety during pregnancy and negative emotionality in early infancy was significant in infants carrying one or more copies of the short allele but not in those homozygous for the long allele. The 5-HTTLPR short allele might increase vulnerability to adverse environmental influences as early as the fetal period.

A genome-wide association study identifies an osteoarthritis susceptibility locus on chromosome 7q22.

OBJECTIVE: To identify novel genes involved in osteoarthritis (OA), by means of a genome-wide association study. METHODS: We tested 500,510 single-nucleotide polymorphisms (SNPs) in 1,341 Dutch Caucasian OA cases and 3,496 Dutch Caucasian controls. SNPs associated with at least 2 OA phenotypes were analyzed in 14,938 OA cases and approximately 39,000 controls. Meta-analyses were performed using the program Comprehensive Meta-analysis, with P values <1 x 10(-7) considered genome-wide significant. RESULTS: The C allele of rs3815148 on chromosome 7q22 (minor allele frequency 23%; intron 12 of the COG5 gene) was associated with a 1.14-fold increased risk (95% confidence interval 1.09-1.19) of knee and/or hand OA (P = 8 x 10(-8)) and also with a 30% increased risk of knee OA progression (95% confidence interval 1.03-1.64) (P = 0.03). This SNP is in almost complete linkage disequilibrium with rs3757713 (68 kb upstream of GPR22), which is associated with GPR22 expression levels in lymphoblast cell lines (P = 4 x 10(-12)). Immunohistochemistry experiments revealed that G protein-coupled receptor protein 22 (GPR22) was absent in normal mouse articular cartilage or synovium. However, GPR22-positive chondrocytes were found in the upper layers of the articular cartilage of mouse knee joints that were challenged with in vivo papain treatment or methylated bovine serum albumin treatment. GPR22-positive chondrocyte-like cells were also found in osteophytes in instability-induced OA. CONCLUSION: Our findings identify a novel common variant on chromosome 7q22 that influences susceptibility to prevalence and progression of OA. Since the GPR22 gene encodes a G protein-coupled receptor, this is potentially an interesting therapeutic target.